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Abstract. We present preliminary results of a long-term 
radial- velocity search for companions to nearby M dwarfs, 
started in September 95. The observed sample is volume- 
limited, and defined by the 127 northern (S > —16°) M dwarfs 
listed in the Gliese and Jahreiss (CNS3) catalogue with d<9pc 
and V<15. Observations are obtained with the ELODIE spec- 
trograph on the 1.93-m telescope of the Observatoire de Haute- 
Provence. The typical accuracy ranges between 10 ms _1 (the 
instrumental stability limit) for the brighter stars and 70 ms 
at our limiting magnitude. We complement the ELODIE veloc- 
ities with older measurements extracted from the CORAVEL 
database to extend our time base, albeit obviously with lower 
precision. Simultaneously, we perform adaptive optics imag- 
ing at CFHT and ESO to look for close (a >0. 05-0.1") visual 
companions in a larger volume-limited sample. For stellar com- 
panions the two techniques together cover the full separation 
range, to beyond the limiting distance of the sample. We will 
therefore eventually obtain a statistically meaningful inventory 
of the stellar multiplicity of nearby M-dwarf systems. We also 
have useful sensitivity to giant planets, as illustrated by our 
recent detection of a planetary companion to Gl 876. 

After 2.5 years, we have discovered 12 previously unknown 
components in this 127 stars sample, plus a companion to an 
additional star beyond its southern declination limit. 7 of these 
are actually beyond the 9 pc limit, as they belong to systems 
included in the sample on the basis of CNS3 photometric paral- 
laxes which were biased-down by the unrecognized companion. 
The remaining 5 companions are true additions to the 9 pc in- 
ventory. More are certainly forthcoming, given our present se- 
lection bias towards short periods and relatively massive com- 
panions. 

We have derived orbital elements for 7 of the new systems, 
as well as for some known binaries with previously undeter- 
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mined orbits. One system, G 203-47, associates an M3.5V star 
with a white dwarf in a rather tight orbit (ai sini = 15 Rq) 
and represents a Post-Common-Envelope system. Some of the 
new M-dwarf binaries will over the next few years provide very 
precise mass determinations, and will thus better constrain the 
still poorly determined lower main-sequence mass-luminosity 
relation. The first such results are now being obtained, with 
some preliminary accuracies that range between 2% at 0.4- 
0.6 Mq and 10% at 0.1 Mq. We have also discovered the third 
known detached M-dwarf eclipsing binary, and determined its 
masses with 0.4% accuracy. 

Key words: stars: binaries - stars: low mass, brown dwarfs - 
planets: giant planets - techniques: radial velocity - techniques: 
adaptive optics 



1. Introduction 

Binarity is a key observational parameter for many astrophys- 
ical questions. One particularly important issue is whether low 
mass stars and brown dwarfs significantly contribute to the 
dark matter in the galactic disk; the present uncertainties on 
the local mass function for very low mass stars are large enough 
for this question to be still unsettled (Tinney 1993), although 
evidence increasingly points towards a negative answer. The 
DENIS and 2MASS near infrared surveys will soon provide 
photometric luminosity functions with very good statistical 
precision, down to the brown-dwarf domain (Delfosse et al. 
1997, 1998a; Kirkpatrick et al. 1997). The corrections for un- 
resolved binaries and the mass-luminosity relation will then be 
the main uncertainty sources for the local stellar mass func- 
tion. Below M~0.3 Mq, the mass- luminosity relation (Henry 
& McCarthy 1993) is only determined by few observational 
data points, most of which are of low accuracy. Accurate orbital 
elements are thus clearly needed for more low mass binaries, to 
provide additional direct stellar mass determinations. A bet- 
ter knowledge of the multiplicity statistics of low mass stars is 
also needed to account for the effect of unresolved systems on 
the mass and luminosity functions: at present, different plausi- 
ble assumptions on stellar multiplicity (Kroupa 1995; Reid & 
Gizis 1997) lead to very different luminosity functions at very 
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low masses. The stellar multiplicity statistics also contains im- 
portant information on both star formation processes and dy- 
namical evolution of stellar systems, as reviewed for instance 
by Duquennoy & Mayor (1991). 

While the binarity statistics are now quite well determined 
for G (Duquennoy & Mayor 1991) and K dwarfs (Halbwachs, 
Mayor and Udry 1998), this is not yet the case for the fainter 
M dwarfs. The only well defined M-dwarf sample which has es- 
sentially complete multiplicity information is the (small) sam- 
ple of M dwarfs within 5.2 pc (Henry & McCarthy, 1990; Lein- 
ert et al., 1997). A number of programmes have searched for 
M-dwarf companions beyond this distance, using different tech- 
niques: radial-velocity monitoring (Bopp & Meredith 1986; 
Young et al. 1987; Upgren & Caruso 1988; Marcy & Benitz 
1989; Tokovinin 1992; Upgren & Harlow 1996), deep infrared 
imaging (Skrutskie et al. 1989; Nakajima et al. 1994; Simons 
et al. 1996), astrometry and speckle or adaptive optics imag- 
ing (Henry & McCarthy 1992; Mariotti et al. 1992). Taken 
together, these programmes are however not sensitive to all 
binary separations for any statistically well defined sample. 
Reid & Gizis (1997) for instance compiled multiplicity infor- 
mation for a northern 8-pc sample, but had to conclude that, 
in spite of Henry & McCarthy's (1992) extensive speckle work, 
the information they could gather remained incomplete. The 
binary fraction derived from these data range from 26% (Lein- 
ert et al., 1997) to 42% (Fisher and Marcy, 1992), through 33% 
(Reid and Gizis, 1997). This seems to points towards a smaller 
fraction of multiple stars than the 57% found by Duquennoy 
& Mayor (1991) amongst G dwarfs, a result which would be 
an important input to multiple star formation theories. All of 
these estimates however suffer from either small number statis- 
tics or large incompleteness corrections, which are quite uncer- 
tain, as they sensitively depend on the (unknown) underlying 
multiplicity distribution. 

Since September 1995, we have therefore been searching 
volume-limited samples of nearby M dwarfs for companions, 
combining three observing techniques which together ensure a 
good sensitivity at all separations for stars and brown dwarfs, 
and some useful sensitivity for giant planets: 

— highly accurate (10-70 ms" 1 ) radial velocities are mea- 
sured with the ELODIE spectrograph at the Observa- 
toire de Haute- Provence (Baranne et al. 1996), which was 
also used in the discovery of the first extra-solar planet 
(Mayor and Queloz 1995); they are combined with older 
CORAVEL data to extend our time base at lower accu- 
racy; 

— high angular resolution near infrared images are obtained 
with the ESO (ADONIS) and CFHT (PUE'O) adaptive 
optics systems, to directly detect relatively massive com- 
panions at small separations; 

— deep adaptive optics coronographic infrared images (Beuzit 
et al. 1997) are obtained with ADONIS at ESO, since 
April 1997; they are sensitive to lower mass companions 
at slightly larger separations, and will be discussed in a 
forthcoming paper. 

In this first paper, we present the radial- velocity pro- 
gramme, discuss its observing and data processing techniques 
in some detail, and present its first results together with some 
complementary adaptive optics data obtained with PUE'O. 
The observed sample and its completeness are discussed in 
Sect. 2 together with the observing strategy, while Sect. 3 
discusses the radial-velocity and adaptive optics observations. 



Section 4 presents the results after 2.5 years, namely the dis- 
covery of 13 new M dwarfs and orbital elements for 7 of the 
new systems. Section 5 discusses the implications of the new 
systems for the solar neighbourhood population. 

2. Sample selection and observing strategy 

2.1. The sample 

The observed radial-velocity sample is presented in detail by 
Delfosse et al. (1998b), who discuss the rotational properties 
of the same stars. Briefly, it contains all M dwarfs listed in 
the third edition of the nearby star catalogue (hereafter CNS3; 
preliminary version, Gliese & Jahreiss, 1991) with a distance 
closer than 9 pc and a B1950.0 declination above —16°. 136 
stars fulfill these criteria among which 7 have to be rejected 
because they are fainter than V=15 (the approximate sensitiv- 
ity limit of the instrument we use), and another 2 because they 
are close companions to much brighter G dwarfs from which 
they cannot be separated by the spectrograph input fiber. At 
9 pc the V=15 limit corresponds to an M6 dwarf, and at 5 pc 
to an M6.5 dwarf. 

We have decided to observe a volume-limited sample be- 
cause we are interested in a fair sampling of the local dwarf 
population, and want to derive unbiased statistical informa- 
tion. The quality of the distances listed in the CNS3 is how- 
ever uneven, so this goal could only be partly reached. The new 
(4 th ) edition of the Yale General Catalogue of Stellar Trigono- 
metric Parallaxes (Van Altena et al. 1995) and the HIPPAR- 
COS catalogue (ESA, 1997) together show that 11 sample stars 
are actually outside the nominal 9-pc sphere, while two omit- 
ted stars are actually within it. The 11 rejected stars have a 
disproportionately large fraction of nearly equal-mass multi- 
ple systems. Most of them were included in the initial sample 
on the basis of CNS3 photometric parallaxes. Their inclusion 
therefore largely reflects the 2 3//2 volume bias incurred in using 
photometric parallaxes for unrecognized equal-mass binaries. 

2.2. Observing strategy 

Since no single observing technique is sensitive to companions 
over all separations, we have chosen to combine imaging and 
radial- velocity techniques, to maximize our sensitivity range. 

As discussed below, the standard errors of the radial- 
velocity measurements range between 10 and 70 ms -1 , de- 
pending on apparent magnitude (and hence largely on spectral 
type, for this volume- limited sample). A 400 ms -1 velocity am- 
plitude is thus always easily detectable, even for the faintest 
stars. A 50 ms -1 velocity amplitude is easily detectable for 
all stars brighter than approximately V=10. Figure [l] shows 
how these conservative minimum velocity amplitudes trans- 
late into minimum detectable companion masses as a function 
of orbital period, for two representative distances and primary 
masses. For a one year orbital period, a 400 ms" 1 velocity am- 
plitude corresponds to a 5 Jupiter mass (Mj) body orbiting a 
0.15 Mq M4V-M5V primary, and 50 ms" 1 corresponds to a 
1.5 Mj body orbiting a 0.4 M Q M2V primary. Even though 
this programme was not chiefly devised to search for planets, 
it therefore has good sensitivity to them, as illustrated by our 
recent detection of a ~2 Mj companion on a 2 month orbit 
around Gl 876 (Delfosse et al. 1998d; also Marcy et al. 1998). 
We will, a fortiori, be sensitive to all stellar and substellar com- 
panions with periods up to a few dozen years. The maximum 
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Fig. 1. Companion detectability in the log(P) versus log(M) plane for our radial-velocity and adaptive optics programmes. The 4 diagrams 
correspond to two representative primary masses of 0.4 Mq (~ M2V, Mjf = 6.3, Baraffe et al. 1998) and O.15M0 (~ M5V, Mk = 8.3, 
Baraffe et al. 1998), at two representative distances of 5 and 9 pc. Typical radial-velocity accuracies are respectively 12 and 70 ms" 1 
for 0.4 Mq and 0.15 Mo stars. The radial- velocity detectability limits correspond to minimum amplitudes of 4cr, or respectively 50 and 
300 ms _1 , and a maximum period of 25 years. Some longer period binaries will be detected, but the observing programme is unlikely to 
last long enough to determine their orbital elements. The solar metallicity models of the E.N.S Lyon group (Allard et al. 1996; Baraffe et 
al. 1998; Allard, Baraffe and Chabrier, private communication) for ages of 1 and 5 Gyr (dashed and solid lines, respectively) are used to 
transform the maximum detectable luminosity contrasts for PUE'O into approximate minimum companion masses. A few representative 
companions of M dwarfs have been added for illustration. Gl 623B and GJ 1245C are two of the least massive objects with dynamical 
mass determinations (Henry and McCarthy 1993). Gl 268. 3B and GJ 2069Ab were discovered in the present programme. The brown dwarf 
Gl 229B (Nakajima et al. 1995) has a most probable mass of 0.05 Mq (Allard et al. 1996). Gl 876b is the first planet discovered around an 
M dwarf (Delfosse et al. 1998d). 



period at which we can actually detect stellar companions will 
be set by the duration of the observing programme rather than 
by the velocity amplitude (Fig. |lj), as it is difficult to detect 
variability on timescales longer than about twice the time-span 
of the data. At 9 pc a limit of P~20 years translates into a min- 
imum angular separation of ~0.4" for a 0.15 Mq primary star 
and a companion of negligible mass (Fig. this separation 
scales as (Mi + M 2 ) 1/3 ). 

The adaptive optics programme (hereafter AO) (Mariotti 
et al. 1992, Perrier et al. 1998) easily resolves equal-mass bi- 
naries down to separations as small as 0.08" (0.7 AU at 9 pc, 
corresponding to a ~1 year period for an M-dwarf pair). It has 
a dynamic range in the K band (Ao = 2.2 /im) of 6.5 mag at a 
0.25" separation, 7.5 mag at 0.5" and 8.5 mag at 1". Since the 
K band absolute magnitudes of M0 and M8 dwarfs only dif- 
fer by ~5 mag (Leggett 1992), the AO observations can detect 



all main-sequence companions that the radial-velocity mon- 
itoring would miss (Fig. |l]). For brown dwarfs on the other 
hand, we do have a small sensivity gap at intermediate separa- 
tions: radial-velocity monitoring can for instance only detect a 
0.05 M Gl 229B-like brown dwarf orbiting a 0.4 M M2V pri- 
mary out to 5.4 AU (0.6" at 9 pc), while AO imaging can only 
detect it beyond 1". This gap disappears for later M-dwarf 
primaries or more nearby systems (Fig. |l|), and the PUE'O 
adaptive optics system for instance easily detects Gl 229B it- 
self, without a coronograph. Beyond 2" (18 AU at 9 pc), the 
coronographic mode of the ESO adaptive optics system (Beuzit 
et al. 1997) can handle an extremely large contrast between 
the primary and a faint secondary (Am~12.5 mag at 2"). We 
are thus again sensitive to all brown dwarf companions of M 
dwarfs, though not to planets. 

Within the [0.4 AU, 4AU] separation range (at 9 pc), the 
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sensitivities of high precision radial velocities and adaptive op- 
tics overlap. In this range, the combination of the two tech- 
niques offers potential access to very accurate masses, as ex- 
emplified by Gl 570B (M1.5V+M3V) for which Forveille et al. 
(1998) obtain masses accurate to ~2%. This should dramati- 
cally improve the mass- luminosity relation for the bottom of 
the main sequence. 

3. Observations 

3.1. ELODIE radial-velocity data 

3.1.1. Instrumental setup 

Most new radial-velocity observations were obtained at the 
Observatoire de Haute-Provence with the ELODIE spectro- 
graph (Baranne et al. 1996) on the 1.93-m telescope, between 
September 1995 and May 1998. This fixed configuration dual- 
fiber-fed echelle spectrograph covers in a single exposure the 
390-680 nm spectral range, at an average resolving power of 
42000. An elaborate on-line processing is integrated with the 
spectrograph control software, and automatically produces op- 
timally extracted and wavelength calibrated spectra, with al- 
gorithms described in Baranne et al. (1996). All stars in this 
programme are now observed with a thorium lamp illuminating 
the monitoring fiber, as needed for the best (~ 10 ms -1 ) radial 
velocity precision (Baranne et al. 1996). Before mid-1997 the 
fainter (V c 13 mag) stars were instead observed with this fiber 
illuminated by the sky, allowing subtraction of the diffused so- 
lar light whose lines could otherwise bias the velocity profile. 
The use of the new M4V correlation template (discussed below) 
has sufficiently reduced the sensitivity of the velocity measure- 
ments to sky emission that the standard high precision setup 
can now be used for all stars in this programme. These early 
data for faint stars have increased random zero-point errors, 
up to ~ 100 ms~ . 

A few measurements were also obtained with the 
CORALIE spectrograph on the newly commissioned 1.2-m 
Swiss telescope at La Silla Observatory (Chile). CORALIE is 
an improved copy of ELODIE and has very similar character- 
istics, with the exception of a substantially improved intrinsic 
stability and higher spectral resolution (R = 50000). 

3.1.2. Data processing 

The extracted spectra are analysed for velocity by digital cross- 
correlation with a one-bit (0/1) template. This processing is 
standard for ELODIE spectra (Queloz 1995a, 1995b), but the 
default K0III mask provided in the ELODIE reduction soft- 
ware is a relatively poor match to the much redder spectra of 
the programme stars. We have constructed a better adapted 
correlation template, applying the method of Baranne, Mayor 
and Poncet (1979) to a high signal-to-noise (S/N=70) ELODIE 
spectrum of Gl 699 (Barnard's stars, M4V). The analysis was 
restricted to the [443 nm, 680 nm] spectral range, since bluer 
orders contain very little flux, for both Gl 699 and the pro- 
gramme stars. For each of the 48 ELODIE orders, 100 trial 
templates go (A) were generated by thresholding the rectified 
and normalized Gl 699 spectrum at 100 levels. We selected for 
each order the thresholding level which maximizes the qual- 
ity factor defined by Baranne, Mayor and Poncet (1979). The 
best order templates were then assembled into one global bi- 



nary template, which has over 3100 "transparent" sections and 
a high overall transmission of 20.75%. 

3.1.3. Measurement accuracy 

For ELODIE, the radial-velocity precision can be written 
(Baranne et al. 1996, their Eq. (9)) as: 

£ (V) - C{T Bff ) (1 + 0.47a) , 
e P [V r )- DS/N 3 kms 

where a is the rotationally broadened 1/e half- width of the 
cross-correlation function (CCF), D its relative depth and S/N 
the signal-to-noise ratio in the reference 48 th ELODIE order. 
The C(T e //) constant depends on both the correlation tem- 
plate used and the spectral type. Its values for the K0 and M4 
templates were determined through Monte-Carlo simulations 
for the [MOV, M5.5V] spectral-type range. High signal-to-noise 
ratio spectra of 6 slowly rotating stars spanning the desired 
range of spectral type (selected from Delfosse et al. 1998b, Ta- 
ble hh were degraded by synthetic photon and readout noise. 
The dispersion of the velocities measured from these spectra 
was then determined for a number of degraded signal-to-noise 
ratio values. The validity of the parameterization was verified 
on those synthetic data, and C(T e ff) was determined through 
a least square fit. It turned out to be constant over the M- 
dwarf spectral range. Its respective values for the K0 (Baranne 
et al. 1996) and M4 templates are 0.085 and 0.035 kms" 1 . The 
radial-velocity standard errors for the two templates are pre- 
sented in Fig. ^. The M4 template improves the accuracy by 
a factor of 2 at spectral type MOV, and by a factor of 10 at 
spectral type M4V. At a given signal-to-noise ratio, the ac- 
curacy variations as a function of spectral type mostly reflect 
changes in the relative depth of the CCF. The width of the 
CCF (Table @) is essentially constant for the K0 template but 
significantly varies with spectral type for the M4 template. 
This change is small and doesn't appreciably affect the preci- 
sion of the radial- velocity measurements, but calibration of the 
CCF width in term of rotational velocity is substantially more 
difficult for the M4 template. 

Table 1. Width and depth of the CCF of 6 stars with spectral 
type in the [M0V,M5.5V] range, for the M4 and K0 templates. 
Spectral types are from Reid et al. (1995). 



Spect. Type 


(JM4(kms x ) 


Dm4 


o~ ko (km s 1 ) 




M0 (Gl 424) 


4.71 


.094 


5.00 


.113 


Ml (Gl 15A) 


4.59 


.104 


4.99 


.106 


M2 (Gl 411) 


4.39 


.119 


4.94 


.096 


M3 (Gl 725A) 


4.32 


.140 


4.87 


.082 


M4 (Gl 699) 


4.21 


.159 


4.96 


.062 


M5 (GJ 1002) 


4.26 


.176 







For the slowly rotating stars in our sample, the above 
formula predicts radial-velocity accuracies that range from 
2 ms" 1 for bright MOV (at typical S/N=200) to 70 ms" 1 for 
faint M6V (at typical S/N=3). It is however clear that the low- 
est values are somewhat optimistic, as the (excellent) long term 
stability of ELODIE radial velocities is ~ 12 ms" 1 (Baranne 
et al. 1996). We have therefore quadratically added 12 ms -1 
to the internal errors to account for instrumental stability. A 
fiber optics light scrambler was recently added to ELODIE 
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Fig. 2. Internal standard errors on radial velocity as a function of signal-to-noise ratio, for the two correlation templates (M4 and KO), 
and for stellar spectral types between MO and M5.5. 



to further stabilize the illumination of the spectrograph, and 
has improved the long term radial-velocity stability to below 
10 ms -1 . We have however chosen not to use this device, as 
its 20% light loss would degrade our overall measurement pre- 
cision for stars fainter than V = 12. 

Under some circumstances, the correction of the radial ve- 
locities to the solar system barycentre represents an additional 
limitation on the accuracy. At the latitude of Observatoire de 
Haute-Provence this correction can vary by up to 1.5 ms -1 
per minute. It therefore significantly changes during the long 
exposure times (up to 1 hour) used for this programme, and 
it is important that the effective time of the observation be 
accurately determined. For bright stars a photomultiplicator 
monitors the fiber illumination of the spectrograph to compute 
its exact value. This P.M. is unfortunately too noisy to be used 
for our programme stars, and we have instead to use the mid- 
point of the observation. Under good atmospheric conditions 
this is not a problem, but when they are unstable (variable see- 
ing or passing clouds) the two values can differ by up to a few 
minutes. In the worst case of unstable atmospheric conditions, 
a 1 hour exposure, and maximum variation of the barycentric 
correction, this additional error can reach 15 ms -1 . It is usu- 
ally much smaller. 

Similar, but much larger, errors can affect measurements 
of very short period spectroscopic binaries (e.g. GJ 2069A, 



discussed below). Their velocities can vary by as much as 
300 ms _1 .min _1 (for P = 1 day), and even small errors on the 
effective time of the observation translate then into significant 
equivalent radial- velocity errors. 

Magnetic activity is another important potential limitation 
for high-precision Doppler- velocity measurements of M dwarfs. 
Many of them have strong chromospheric activity, visible as 
H a and Ca + H and K emission lines. The surface temperature 
of active stars is inhomogeneous, and since hotter parts are 
brighter, the velocity measured from the disk-integrated spec- 
trum is biased towards the velocity of any hot spot, or away 
from the velocity of any cold spot. This effect is exacerbated 
when using a mismatched correlation template, such as the KO 
template for M dwarfs: the KO template is a better match to 
the hotter parts of the stellar surface and it therefore further 
increases their relative weight in the measured velocity. 

On Julian day 2449978.5 (18.09.1995), a strong chromo- 
spheric flare in Gl 873 provides an extreme example of activity- 
induced radial-velocity errors, and illustrates the much better 
immunity of the M4 correlation template. This generally active 
star presented on this date an activity well above its usual level. 
The H a emission line showed a strong broad pedestal, charac- 
teristic (e.g. Jones et al. 1996) of active stars observed within 
about an hour of a major flare, and usually interpreted (Eason 
et al. 1992) as arising from mass motion of gas clumps with a 
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large volume rilling factor. With the KO template, the radial 
velocity of Gl 873 for this date has a precision of 150 ms -1 but 
differs from the median Gl 873 value by 1.7 kms -1 . With the 
M4 template on the other hand, the radial velocity precision for 
the same spectrum is improved to 15 ms _1 , and the measured 
value is within only 40 ms _1 of the median velocity. Given 
the extreme character of the GL 873 event, 40 ms _1 can thus 
be taken as an upper bound on the maximum radial-velocity 
error that can result from unusually high magnetic activity, 
when using the M4 mask. 

The examination of the ~50 apparently single stars which 
have at least 3 good measurements (internal precision better 
than 20 ms _1 ) shows that about half of them have apparent 
excess velocity dispersion at the ~ 30 ms^ 1 level. Gl 876 has 
been one of those until we very recently accumulated enough 
measurements to ascribe its velocity variability to a planetary 
companion (Delfosse et al. 1998d; also Marcy et al. 1998). Most 
of these stars however have few measurements, and it is thus 
not yet possible to determine how much of the apparent dis- 
persion is due to intrinsic or instrumental radial-velocity noise, 
and how much is due to yet unrecognized companions, either 
planetary ones in short-period orbits or more massive distant 
ones. The exact accuracy limits of this programme are thus 
still undetermined, but are better than ~ 30 ms _1 . 

Double-lined spectroscopic binaries are in addition affected 
by a phase-dependent systematic error source. The Gaussian 
profile adjustment which is used to measure the radial veloc- 
ity only imperfectly describes the correlation profile of an M 
dwarf, whose baseline shows some low level (~1%) systematic 
undulations. For a single star, or a single-lined binary, this has 
no effect on the measured velocity, or perhaps introduces a 
small constant zero-point offset. For double-lined binaries on 
the other hand, the correlation peak of one star is superposed 
onto a systematically different part of the other star's base- 
line, producing phase-dependent velocity errors. If the velocity 
amplitudes are of the order of the width of the the two corre- 
lation profiles, these are systematic errors and they will bias 
the derived orbital elements, and in particular the amplitudes. 
When the amplitudes are sufficiently larger than the width of 
the profiles on the other hand, these velocity errors become 
random, and the orbital elements then only have increased un- 
certainties but are unbiased. Under the same circumstances, 
and with enough measurements, it is also possible to separately 
determine the correlation profiles of the two components, and 
to correct the velocities for the effect of their baseline wig- 
gles, for instance with the iterative method used by Forveille 
et al. (1998) for Gl 570B. For smaller velocity amplitudes this 
method cannot be used, and one would have instead to rely on 
the correlation profile of a matched spectral template, intro- 
ducing some additional uncertainties. We haven't yet applied 
any such correction except to Gl 570B, and SB2 velocities thus 
have errors of ~ 100 ms _1 or larger (for the fainter component 
in particular). 

3.2. CORAVEL radial-velocity data 

65 bright (V < 12) stars in the sample had previously been 
measured in the course of various observing programmes using 
the older CORAVEL spectrographs (Baranne et al. 1979) on 
the 1-m Swiss telescope (at Observatoire de Haute-Provence, 
France) and the 1.54-m Danish telescope (La Silla Observatory, 
ESO, Chile). The number of measurements varies considerably, 



from 2-3 for stars examined in the course of the HIPPARCOS 
follow-up survey (Udry et al. 1997), to over 100 for a few well 
studied previously known binaries. They have lower precision 
than the ELODIE measurements (at best 300 ms _1 , and often 
1-2 kms -1 for the present faint red stars), but, when available, 
they extend the time base to typically 10 years. These data are 
thus systematically extracted from the CORAVEL database 
for all programme stars. There is a significant radial velocity 
offset between CORAVEL and ELODIE measurements (Udry 
et al. 1998), which must be corrected before they can be com- 
bined. For G-K dwarfs the offset depends upon both spectral 
type and radial velocity. Over the [MOV, M3.5] range where we 
have common measurements the offset is independent of spec- 
tral type, but has a significant uncalibrated star to star scatter 
at the ~ 600 ms -1 level. For well observed binary stars, the 
CORAVEL offset is thus determined as part of the orbital el- 
ement adjustment. 

3.3. Adaptive optics data 

3.3.1. Instrumental setup 

The observations were carried out at the 3.6-meter Canada- 
France-Hawaii Telescope (CFHT) during several observing 
runs from September 1996 to March 1998, using the CFHT 
Adaptive Optics Bonnette (AOB) and two different infrared 
cameras. The AOB, also called PUE'O, is a general pur- 
pose adaptive optics (AO) system. It is mounted at the F/8 
Cassegrain focus, and cameras and other instruments are then 
attached to it (Arsenault et al. 1994, Rigaut et al. 1998). It 
analyses the atmospheric turbulence with a 19 element wave- 
front curvature sensor and corrects for it with a 19 degree of 
freedom bimorph mirror. Modal control and continuous mode 
gain optimization (Gendron & Lena 1994; Rigaut et al. 1994) 
maximize the correction quality for current atmospheric turbu- 
lence and guide star magnitude. For our observations a dichroic 
mirror diverted the visible light to the wavefront sensor while 
a science detector recorded near- infrared light. The data were 
recorded with either MONICA, the Universite de Montreal In- 
frared Camera (Nadeau et al. 1994), or KIR, a new CFHT 
infrared camera developed to take full advantage of the AO 
corrected images produced by PUE'O (Doyon et al. 1998). 

MONICA was used for the commissioning of the AOB dur- 
ing the first semester 1996 and for all science runs until Novem- 
ber 1997. It is a facility instrument based on a NICMOS3 256 x 
256 detector, and was originally designed by the Universite de 
Montreal for the Observatoire du Mont Megantic and CFHT 
F/8 Cassegrain focii. The camera was refitted with new optics 
for use at the F/20 output focus of AOB. It produces a plate 
scale of 0'.'034 per pixel, properly sampling diffraction-limited 
images down to the J band (1.25 /xm). The resulting field size 
is 8.7" x 8.7". 

Since December 1997, MONICA has been replaced on 
PUE'O by KIR, a new imaging camera which records a 16 
time larger field on an HAWAII 1024 x 1024 HgCdTe array. 
KIR has improved optical quality and detector read-out noise, 
and therefore a significantly better detectivity. The KIR plate 
scale is 0"035 per pixel, for a total field size of 36" x 36". 

3.3.2. Observations 

Sources were first examined for binarity with one filter, usually 
H (1.65 nva). Under good to moderate seeing conditions H rep- 
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resents the best compromise between sensitivity, corrected im- 
age quality, and sky brightness. Under worse seeing conditions 
the K filter was used instead to maintain acceptable image 
quality. Sources which saturate the detectors in the minimum 
available integration time through the H (1.65 lira) or K(2.23 
pm) broad-band filter (brighter than K — 7 under typical con- 
ditions) were observed through corresponding narrow-band fil- 
ters, respectively [Fe + ] (1.65 (im) and Hh (2.12 fim). Whenever 
a target appeared double, it was further observed with addi- 
tional filters to determine relative colour indices. Integration 
times per frame typically range between a few tenths of a sec- 
ond and a few seconds. In order to improve the signal-to-noise 
ratio and to average the residual uncorrected atmospheric tur- 
bulence, series of ~4 minute total integration times were accu- 
mulated in a four position mosaic pattern. This observing se- 
quence also allows to correctly determine the sky background 
and to correct for detector cosmetic defects. Wavefront sens- 
ing was performed on the sources themselves, which are always 
bright enough (R < 14) to ensure diffraction-limited images in 
H and K bands under standard Mauna Kea atmospheric con- 
ditions (i.e. seeing up to 1"). For most targets, PSF calibration 
stars were observed under the same conditions to provide input 
to parameter fitting (Section 3.3.3| ) and image deconvolution. 
The atmospheric turbulence and AO correction for a given set 
of observations were further characterized by simultaneously 
recording the wavefront sensor measurements and deformable 
mirror commands. An accurate synthetic PSF can be gener- 
ated a posteriori from these ancillary data, as described by 
Veran et al. (1997). Astrometric calibration fields such as the 
central region of the Trapezium Cluster in the Orion Nebula 
(McCaughrean and Stauffer 1994), were observed to accurately 
determine the actual detector plate scale and position angle 
(PA.) origin. Flat-fields were obtained on the dome and the 
sky for each filter. 



3.3.3. Data reduction 

For each filter, the raw images were median combined to pro- 
duce sky frames which were then subtracted from the raw data. 
Subsequent reduction steps included flat-fielding, correction 
from the bad pixels, and finally shift-and-add combinations 
of the corrected frames into a final image. For resolved binary 
systems, the separation, position angle and magnitude differ- 
ence between the two stars were determined using uv plane 
model fitting in the GILD AS (Grenoble Image and Line Data 
Analysis System) software. With approximate initial values of 
the positions of the two components along with a PSF refer- 
ence image, the fitting procedures gave as output the flux and 
pixel coordinates of the primary and secondary. Application 
of the astrometric calibrations then yields the desired param- 
eters. Images of the newly resolved binaries are presented in 

Fig. |. 



4. New companions 

In this section we comment the new companions individually. A 
general view of their measurements (visual and spectroscopic) 
and status (new discovery, first orbit determination, etc.) is 
summarized in Table ^| If determined the orbital elements are 
given in Table ^ and some interesting visual parameters in 
Table |. 



Table 2. New low mass companions in the solar neighbour- 
hood (the companion of Gl 876 is a planet). Spectral types are 
from Reid et al. (1995), except for GJ 2130B which is taken 
from the CNS3 catalogue. Parallaxes (in milliarcseconds) are 
taken from (a) the HIPPARCOS catalogue, (b) the Yale cata- 
logue (Van Altena et al. 1995) or (c) are photometric parallaxes 
from the CNS3. 



Name 


parallax (mas) 


joint 






spectral type 


LP 476-207 AabB 


91.20±8.56 (a) 


M 


4 


Gl 268.3 AB 


81.05±2.42 (a) 


M 


2.5 


GJ 2069 Aab 


78.05±5.69 (a) 


M 


3.5 


GJ 2069 BC 


78.05±5.69 (a) 


M 


4 


LHS 6158 AabB 


224.0±36.0 (c) 


M 


3.5 


Gl 381 AB 


81.23±2.37 (a) 


M 


2.5 


Gl 487 AabB 


98.14±1.67 (a) 


M 


3 


LHS 2887 AB 


62.2±13.1 (6) 


M 


4 


G 203-047 ab 


137.84±8.95 (a) 


M 


3.5 


GJ 2130 Bab 


161.77±11.29 (a) 


M 


2.5 


Gl 829 ab 


148.29±1.85 (a) 


M 


3.5 


Gl 876 ab 


212.69±2.10 (a) 


M 


4 


Gl 896 Aab 


160.06±2.81 (a) 


M 


3.5 


Gl 896 Bab 


160.06±2.81 (a) 


M 


4.5 



4.1. LP 476-207 AabB 

LP 476-207 is a new triple system. Adaptive optics images 
show a 0.97" separation pair with a K-band magnitude dif- 
ference of 0.9 (Table |), confirming a recent speckle detection 
of this outer component by Henry et al. (1997). The brightest 
component of the visual pair, LP 476-207 A, is a new large 
amplitude double-lined spectroscopic binary, with a period of 
P~12 days. The visual pair was also detected by HIPPARCOS 
(ESA, 1997), which lists a separation of 0.68". There has thus 
been significant orbital motion over the last ~5 years. 

4.2. Gl 268 ab 

Gl 268 is a previously known double-lined binary on a 10.43- 
day orbit (Tomkin and Pettersen 1986) . Our new spectroscopic 
measurements provide substantially improved orbital elements. 

4.3. Gl 268.3 AB 

Gl 268.3 is a new double-lined spectroscopic binary. With re- 
spectively 29 and 10 CORAVEL and ELODIE measurements, 
spanning over 15 years, the radial-velocity orbit is extremely 
well determined. The period is 304.35±0.25 days. K-band 
adaptive optics images partly resolve this system into a ~0.1" 
binary with a small magnitude difference (Table |B|). Within 
one or two years Gl 268.3 will thus provide two precise mass 
determinations at a spectral type of M2.5V. 

44. GJ 2069 AabBC 

Previously known as a wide (~ 12") visual binary, the GJ 2069 
system is actually quadruple. Adaptive optics images resolve 
GJ 2069B into a 0.36" binary with a K-band magnitude dif- 
ference of 0.45 (Table |^). The 5 radial- velocity measurements 
obtained over ~850 days show a linear drift of ~ 600 ms _1 . 
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Fig. 3. Radial-velocity orbits. Triangles correspond to radial velocities of the primary star and squares to those of the secondary star. 
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The period is thus probably long and the small visual separa- 
tion may be due to projection. 

GJ 2069A is a new short-period (P = 2.8 days) detached 
double-lined eclipsing binary, only the third one known with 
substantially sub-solar masses. Delfosse et al (1998c) discuss in 
detail the very accurate radial-velocity orbit and a preliminary 
light curve. They derive individual masses (0.430 and 0.396 
Mq) with 0.4% accuracy. Both components are clearly sub- 
luminous for their masses, and their spectral types are later 
than expected, compared with both solar metallicity stellar 
models and previous stellar mass determinations. This indi- 
cates a metal-rich composition. 

4.5. LHS 6158 (G 041-014) AabB 

Recently shown to be a short-period double-line spectroscopic 
binary (Reid & Gizis 1997), G 41-14 is actually a triple system. 
We have determined the orbital elements (P — 7.6 days) of 
the close spectroscopic pair, and discovered a third component 
in AO images, at a separation of 0.62" and with a K-band 
magnitude difference of 0.5 (Table ^). The probable period 
is ~10 years and a mass determination is thus a mid-term 
prospect only. 

4.6. Gl 381 AB 

This new long-period double-line spectroscopic binary is also 
well resolved into a 0.18" pair (AK=0.95 mag) in AO im- 
ages (Table |). The period is relatively well determined at 
~2845 days, but the CORAVEL measurements do not sepa- 
rate the two components well and the ELODIE data only cover 
a small fraction of the orbit. All other orbital elements are thus 
still uncertain. ELODIE will easily separate the two stars at 
periastron. This spectroscopic + visual pair will then provide 
accurate mass measurements, within about 5 years. 

4.7. Gl 487 AabB 

Gl 487 is a new triple-line spectroscopic system. The short 
period (for Gl 487Aab) is 54.07 days and the long period is 
approximately 3000 days. The outer system (Gl 487 AB) is 
also resolved in AO images with a separation of 0.23" and 
AK=0.7 mag (Table ^). This system has the potential to pro- 
vide three accurate mass determinations, but its complete anal- 



ysis is deferred to a forthcoming paper (Delfosse et al. in prepa- 
ration) dealing with triple systems. 

4.8. LHS 2887 (G165-061) AB 

The probable period of this new double-line spectroscopic bi- 
nary is a few years (1500 days for a very preliminary orbit). It 
is also well resolved in AO images, with a separation of 0.58" 
(Table |) and will eventually provide accurate masses. 

4.9. Gl 644 

Our spectroscopic measurements confirm that the speckle bi- 
nary Gl 644 (Tokovinin et Ismailov 1988) is a triple system, as 
previously shown by Eggen (1978) and Pettersen et al. (1984). 
We have decoupled the two orbits, and obtained preliminary 
values for the masses of the three stars. Their accuracy would 
however benefit from additional measurements and their dis- 
cussion is thus deferred to a forthcoming paper (Delfosse et al., 
in preparation). 

4.10. G 203-047 ab 

This M3.5V star was previously noted as a single-line spec- 
troscopic binary by Reid & Gizis (1997). We have obtained 
orbital elements, with a 15-day period, a large velocity semi- 
amplitude of 50 kms -1 , and a small and only marginally sig- 
nificant eccentricity of 0.07. It is listed in the HIPPARCOS 
catalogue as a probable short-period astrometric binary with- 
out orbital solution. Now knowing the period, it would be of 
interest to reanalyse the HIPPARCOS data, and extract the 
inclination. The orbital elements result in a large mass func- 
tion, (M 2 sini) 3 /(A/i + A/ 2 ) 2 = 0.2 M , and therefore imply 
a quite massive secondary. Adopting the Baraffe et al (1998) 
theoretical mass- luminosity relation, the mass of the M3.5V 
primary (Mi=8.87, Figueras et al. 1990) is 0.30-0.35 M Q , and 
the measured mass function thus requires that the secondary 
component be more massive than ~O.5M0. This clearly ex- 
cludes a single-main sequence star, which would have a spec- 
tral type earlier than M1V and would dominate the observed 
light. A (very) short-period main-sequence binary is similarly 
excluded, because its brighter member would still need to be at 
least about as massive as the primary, and would be easily visi- 
ble in the spectrum. The companion to G 203-047 must thus be 
a degenerate star. It is most likely a white dwarf, since visible 



10 



Table 3. Measurement summary: COR, ELO and PUE'O are the number of measurements obtained with CORAVEL, ELODIE 
and PUE'O, respectively. Gl 644 and Gl 866 are not included here, since their discussion is deferred to a forthcoming paper 
(Delfosse et al. in preparation). To date (September 1998), this programme has identified 13 new stellar companions and 1 
planet orbiting M dwarfs of the solar neighbourhood. Orbital elements have been determined for six of them, as well as for 3 
previously known binaries with undetermined orbits (LHS 6158 Aab and G 203-047 ab, Reid and Gizis (1997); GJ 1230 Aab, 
Gizis and Reid (1996)). 



Name 


COR 


ELO 


PUE'O 


discovered by 


first orbit 


orbit 


In the 9-pc 










this programme 


determination 


improvement 


sample 


LP 476-207 Aab 





16 


3 


yes 


yes 




no 


LP 476-207 AB 





16 


3 


no 


no 


no 


no 


Gl 268 ab 





9 


1 


no 


no 


yes 


yes 


Gl 268.3 AB 


29 


10 


3 


yes 


yes 


- 


no 


GJ 2069 Aab 





18 


1 


yes 


yes 


- 


no 


GJ 2069 BC 





5 


1 


yes 


no 


no 


no 


LHS 6158 Aab 





14 


2 


no 


yes 


- 


yes 


LHS 6158 AB 





14 


2 


yes 


no 


no 


yes 


Gl 381 AB 


27 


4 


2 


yes 


no 


no 


no 


Gl 487 Aab 


24 


18 


5 


yes 


yes 




no 


Gl 487 AB 


24 


18 


5 


yes 


no 


no 


no 


LHS 2887 AB 





5 


2 


yes 


no 


no 


no 


G 203-047 ab 





12 


1 


no 


yes 




yes 


GJ 2130 Bab 





1 





yes 


no 


no 


yes 


GJ 1230 Aab 





11 





no 


yes 




yes 


Gl 829 ab 


67 


11 


1 


yes 


yes 




yes 


Gl 876 ab 





32 


1 


yes 


yes 




yes 


Gl 896 Aab 





9 


1 


yes 


no 




yes 


GL 896 Bab 





6 


1 


yes 


no 




yes 



Table 4. Orbital elements and Msini for the newly determined radial- velocity orbits. The inclination of GJ 2069 Aab is 
accurately known and we thus list the actual masses instead of Msini. G203-47 and Gl 876 are single-lined binaries (all others 
are double-lined), so we list fi(M) as well as estimated M2 sin i for assumed primary masses of, respectively, 0.35 and 0.30 Mq. 



Name 


P 


To 


e 


UJ 


Ki 


K 2 


V 


mass determination 




(days) 


(Julian day) 






(kins" 1 ) 


(kms- 1 ) 


(kms -1 ) 


(M 


a) 


Double-lined binaries: Mi x sin i and M2 x 


sin i 










Mi sin^ i 


M2 sin J i 


LP 476-207 Aab 


11.9623 


49799.47 


.323 


212.0 


9.96 


17.57 


17.26 


.014 


.008 




±.0005 


±0.04 


±.006 


±0.6 


±0.03 


±0.07 


±0.03 


±.0001 


±0.0001 


Gl 268 ab 


10.4265 


50149.902 


.321 


212.1 


34.81 


40.86 


41.83 


.215 


.183 




±.00002 


±0.008 


±.001 


±0.3 


±0.04 


±0.06 


±0.03 


±.001 


±.001 


Gl 268.3 AB 


304.35 


48826.0 


.399 


273.8 


12.47 


18.6 


-7.51 


.435 


.292 




±.25 


±1.5 


±.008 


±0.9 


±0.08 


±0.2 


±0.05 


±.010 


±.006 


LHS 6158 Aab 


7.5555 


50471.2 


.014 


7.0 


30.15 


36.79 


-7.57 


.129 


.106 




±.0002 


±0.2 


±.002 


±9.3 


±0.05 


±0.09 


±0.04 


±.001 


±.001 


Gl 487 Aab 


54.075 


50506.2 


.081 


137.0 


22.0 


22.5 


-10.1 


.247 


.242 




±.006 


±0.5 


±.005 


±3.0 


±0.5 


±0.5 


±0.4 


±.014 


±.013 


GJ 1230 Aab 


5.06880 


50643.7 


.009 


230.0 


46.9 


49.0 


-11.88 


.237 


.226 




±.00005 


±0.2 


±.001 


±10.0 


±0.1 


±0.1 


±0.05 


±.001 


±.001 


Gl 829 ab 


53.221 


48980.2 


.374 


300.0 


18.7 


18.7 


-25.23 


.114 


.114 




±.004 


±0.2 


±.004 


±1.0 


±0.1 


±0.1 


±0.06 


±.001 


±.001 


Single-lined binaries: /1 (M) and an estimate of M2 sin 










MM) 


M 2 sin a i 


G203-47 Aab 


14.7136 


50500.8 


.068 


175.0 


50.6 




-9.7 


0.2 


~ 0.5 




±.0005 


±0.1 


±.004 


±3.0 


±0.2 




±.2 






Gl 876 ab 


61.1 


50661.7 


.33 


5.0 


0.247 




-1.901 


8.10 s 


~ 0.002 




±.2 


±1.5 


±.02 


±5.0 


±0.006 




±0.005 






Double-lined eclipsing binary: complete determination of Mi and M2 


Mi 


M 2 


GJ 2069 Aab 


2.771472 


50207.8128 


.0 




68.03 


73.06 


4.36 


.430 


.396 




±.000004 


±0.0009 


±.003 




±0.09 


±0.09 


±0.04 


±.001 


±.001 
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photometry of G 203-47 (Figueras et al. 1990) shows a ~0.4 U- 
B excess over the colour of stars with the same R-I colour or 
spectral type (Leggett 1992). UV spectroscopy would easily as- 
certain the exact characteristics of the white dwarf. The small 
semi- major axis of the present orbit (aisini = 14.5 R Q , or 
0.05 AU) implies that G 203-47a must have been, at previous 
stages, in a contact configuration with the AGB progenitor 
of its white dwarf companion. It is therefore a very nearby 
member of the pre-cataclysmic variable family (e.g. Ritter & 
Kolb 1998), but the timescales for its evolution into a CV is 
extremely long. It should have accreted some nucleosynthesis 
products dredged up to the surface of its previous AGB com- 
panion, and its composition may thus be peculiar. 



4.11. GJ 2130 ABab 

At 5 = —32°, the little studied GJ 2130 system is far below the 
declination limit of the main sample and was observed as part 
of a possible southern extension of the programme. GJ 2130B, 
the fainter component of this wide (25" separation) visual bi- 
nary turns out to be a double-lined spectroscopic binary, and 
the system is thus triple. To date we don't have enough radial- 
velocity data to attempt an orbital solution. 

4.12. GJ 1230 AabB 

GJ 1230A, the brighter component of the wide (5" separation) 
GJ 1230AB visual binary is a double- lined spectroscopic binary 
(Gizis and Reid, 1996), for which we determine first orbital 
elements. The period is 5.1 days. 

4.13. Gl 829 ab 

Gl 829 was mentioned as a possible double-lined spectroscopic 
binary by Marcy et al. (1987). It is clearly seen as such by 
both CORAVEL (67 measurements) and ELODIE (11 mea- 
surements) and the orbital elements are well constrained. The 
orbital period is 53.2 days and the mass ratio is very close to 
1. 

4.14. Gl 866 

We confirm that the astrometric and speckle binary Gl 866 
(Leinert et al. 1986) is a triple system, as long suspected from 
the mass excess of one of the speckle components (Leinert et 
al. 1990). It is seen as a triple-lined spectroscopic system in the 
ELODIE spectra. The extensive speckle coverage of its outer 
orbit (Leinert et al. 1990) makes it a prime candidate for an 
accurate mass determination. As for Gl 644, the two decoupled 
orbits and the resulting masses will be discussed in a forthcom- 
ing paper (Delfosse et al., in preparation). 

4.15. Gl 876 ab 

High precision radial-velocity observations of the nearby M4 
dwarf Gl 876 with the Observatoire de Haute-Provence 1.93-m 
telescope and the new 1.20-m Swiss telescope at La Silla indi- 
cate the presence of a Jovian mass companion to this star. The 
orbital fit to the data gives a period of 60.96 days, a velocity 
amplitude of 248 ms^ 1 and an eccentricity of 0.34. Assuming 
that Gl 876 has a mass of 0.3 M Q , the mass function implies a 
mass for the companion of 2/sinj Jupiter masses. Delfosse et 



al. (1998d) and Marcy et al. (1998) both discuss this interesting 
system in detail. 

4. 16. Gl 896 AabBab 

The Gl 896 system, well known as a wide binary, is actually 
quadruple: both members of the visual pair are new single- 
lined spectroscopic binaries. Their probable periods are of the 
order of a few years but still undetermined. We have observed 
both of them with adaptive optics in December 1997 and not 
resolved either. The luminosity contrast of the companions is 
thus probably large. 



Table 5. Projected separation and luminosity contrast at K 
band, for binaries resolved by adaptive optic imaging. 



Name 


Separation 


AK 

(mag) 


Observing 
date 


LP 476-207 AB 
Gl 268.3 AB 
GJ 2069 BC 
LHS 6158 AB 
Gl 381 AB 
Gl 487 AB 
LHS 2887 AB 


0.97" ±0.01" 
marginally re 
0.36" ±0.01" 
0.62" ±0.01" 
0.18" ±0.02" 
0.23" ±0.01" 
0.58" ±0.01" 


0.9±0.05 
solved ~0.1" 
0.45±0.3 
0.50±0.04 
0.95±0.05 
0.7±0.1 
0.15±0.03 


Jan 1997 
May 1997 
May 1997 
Jan 1997 
Aug 1997 
Feb 1997 
May 1997 



5. Discussion 

After 2.5 years, we have found 12 new companions to the 
127 M dwarfs in our initial sample, and over the same pe- 
riod 6 additional companions were discovered by others (we 
also independently found most of those). While showing that 
the multiplicity information was previously fairly incomplete, 
this should not be taken as the final number of companions 
for this sample. Many stars still have few radial-velocity mea- 
surements (the median is 5 but a number of stars only have 
2 measurements), which span an interval of at most 2.5 years. 
Our detectivity is thus very significantly biased towards nearly 
equal-mass binaries (detectable as double-lined binaries in a 
single measurement) or/and short periods. A significant num- 
ber of lower mass companions in wider orbits certainly remains 
to be found in this sample. 

The discovery of this large number of new multiple systems 
has two opposite effects on estimates of the stellar density in 
the solar neighbourhood. On the one hand, the new compo- 
nents within 9 pc increase the previously underestimated local 
density, typically with lower mass stars as they are necessar- 
ily fainter than their primaries. On the other hand, a number 
of previously unknown M-dwarf binaries are actually beyond 
the 9-pc distance limit of the sample, in which they were ini- 
tially included on the basis of an underestimated photometric 
distance, reflecting the well known ~2 3 ^ 2 volume bias in photo- 
metric parallaxes of unrecognized binaries. This second correc- 
tion decreases the local density, more or less uniformly for all 
masses. These two effects have been abundantly discussed in 
the context of photometric luminosity functions, as a possible 
explanation of their differences from the solar neighborhood 
luminosity function (e.g. Kroupa 1995, and Reid & Gizis 1997, 
for two contrasted views). It is perhaps not always recognized 
that, beyond ~5 pc, the solar neighborhood luminosity func- 
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Fig. 4. Adaptive optics images for 5 of the new binaries. 



tion still has a significant photometric component, and that it 
is thus also affected at some level by the same two effects. 

Taking together newly published trigonometric parallaxes 
(ESA, 1997; Van Altena at al. 1995) and corrected photomet- 
ric distances to the new binaries, 13-M-dwarf systems listed 
within 9 pc in the preliminary version of the CNS3 are actu- 
ally beyond this distance, while only one system (Gl 203) enters 
this volume. The better observed 5.2-pc sample, on the other 
hand, is essentially unaffected, with only the G041-014 sys- 
tem removed (~4.5 pc — ► ~8 pc) from it. The previously men- 
tioned (e.g. Henry et al. 1994) incompleteness of the nearby 
M-dwarf system sample beyond 5 pc is thus made more signif- 
icant. Systems with bright secondaries are rapidly recognized 
as double-lined binaries, and their inventory should thus be 
essentially complete. The remaining 7 photometric parallaxes 
in this revised 9-pc sample can therefore be considered as re- 
liable. The 9-pc northern sample of systems should thus now 
only significantly change through inclusion of presently missing 
systems. 

The new multiple systems, on the other hand, add 6 stars to 
the 9-pc sample of stars. Given the remaining selection biases 
in our binary search, we feel that derivations of the binarity 
statistics and stellar luminosity function within 9 pc would still 
be premature. 

Our sensitivity to lower companion masses and longer peri- 
ods is however quickly improving, so that this information will 



soon become available. It will then be possible to more reliably 
estimate the correction for unresolved binaries in photometric 
determinations of the luminosity function, and hopefully settle 
the long-standing controversy (e.g. Kroupa 1995 and Reid & 
Gizis 1997) on its true importance. 
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